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Precision of Infiaton Potential Reconstruction from CMB 
Using the General Slow-Roll Approximation 
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Through a principal component analysis, we study how accurately CMB observables can 
constrain infiaton potentials in a model independent manner. We apply the general slow-roll 
approximation in our analysis where we allow, in contrast to the standard slow-roll approximation, 
the possibility of variations in V"(<j>) and take into account the fact that horizon crossing is not an 
instantaneous event. Our analysis provides a set of modes to be used in fitting observables. We 
find that of order five of these modes will be constrained by future observations, so a fully general 
data analysis package could use the amplitudes of just a handful of modes as free parameters and 
retain all relevant information in the data. 

PACS: 98.80.Cq 



I. INTRODUCTION 

The origin of cosmic density perturbations is one of 
the most notable predictions of the inflationary scenario. 
The observational data on their power spectrum from the 
cosmic microwave background (CMB) and large scale 
structure provide us with one of few windows through 
which we can probe the physics governing the early Uni- 
verse [j], 0, 0, Q| ■ In particular it is interesting to ask if 
the CMB can directly constrain the potential energy of 
the field driving inflation, the infiaton, and if so, which 
features of the potential will be most accurately mea- 
sured. Extracting the physics of inflation from CMB 
data alone is not trivial because different inflation models 
can produce identical CMB spectra and moreover CMB 
observables can cover only a tiny portion of an infiaton 
potential. However, we can still try to constrain the cur- 
vature, for example, of the infiaton potential during the 
time when cosmologically interesting scales were exit- 
ing the horizon, and the sign of the curvature alone can 
already distinguish between several classes of inflation 
models. Because we are still not certain of the func- 
tional form of the infiaton potential, it would be useful 
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to work in a model-independent manner to avoid bias- 
ing our estimates of the properties of infiaton potentials. 
Further, we are motivated to ask what other features of 
the potential are constrained by the CMB, again in a 
model-independent way? 

There exist many works which reconstruct the primor- 
dial power spectrum P(k) from the observable CMB data 
Ci (see, for example, |1 B H H El E3 for model- 
independent approaches). If we can also measure the 
tensor perturbations and assume the standard slow-roll 
approximation, reconstructing the infiaton potential is 
feasible [12|. Even without the standard slow-roll ap- 
proximation, the conversion from V(k) to infiaton po- 
tcntial parameters is possible analytically for a given 
V(k) 0, 0| . So we could combine these analyses to 
examine the properties of inflation from Ci by two steps: 
Ci — ► V(k) — ► infiaton potential parameters. However, 
degeneracies between the parameters arise from these 
two steps, and it would be more efficient to constrain 
an infiaton potential directly from C/. If we perform 
a numerical analysis such as the Markov Chain Monte 
Carlo method, reconstructing the infiaton potential di- 
rectly from C\ can be done without heavily relying on 
the standard slow-roll approximation p^Lll6|. However, 
even in the case of numerical analysis, how many in- 
fiaton potential parameters do we need in addition to 
cosmological parameters such as the reionization optical 
depth and baryon density? Certainly including more pa- 
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rameters would lead to a better fit to observations, but 
adding as many free parameters as we want would not be 
a practical approach. Our principal component analysis 
will be useful in choosing an optimal set of parameters 
to characterize the inflaton potential. 

In Sec. [n] we describe the general slow-roll approx- 
imation with an emphasis on the difference from the 
standard slow-roll approximation, and discuss a model- 
independent parametrization of the inflaton potential. 
Sec. IIIII outlines a principal component analysis and 
determines the precision that idealized future CMB ex- 
periments can attain in constraining inflaton potential 
parameters. Sec. II VI gives some illustrations of the for- 
malism, followed by the conclusion in Sec. [V] 



II. INFLATON POTENTIAL 
PARAMETRIZATION 

We, in this section, discuss a model independent 
parametrization of the inflaton potential (strictly speak- 
ing, a function of the inflaton potential) in view of the 
general slow-roll approximation to study how precisely 
CMB observables can ultimately constrain the inflaton 
potential. 



A. General Slow- Roll Approximation 

We refer the readers to the Appendix and Ref. 17] for 
the discussion of general slow-roll approximation, but let 
us briefly here point out the motivations for our using the 
general slow-roll approximation rather than the standard 
slow-roll approximation in our analysis. 

One of the possible pitfalls of the standard slow-roll 
approximation is that it presumes the scale invariance of 
V"((j)) (in the rest of the paper, primes denote deriva- 
tives with respect to the argument). For single field in- 
flation with slow-roll, the observations are indeed con- 
sistent with a small amplitude 1 of V". Neither obser- 
vations nor theories, however, require exact scale invari- 
ance. The general slow-roll approximation can lift this 
extra (and unnecessary) condition of the scale invariance 
of V" 17]. This would help us consider the possibilities 
for a wider variety of inflaton potential properties from 
CMB observables. The application of the general slow- 
roll approximation in our principal component analysis 
where the standard slow-roll approximation is not appli- 
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cable is given in the second example of Section IIVI for 
the illustration purpose. 

Another pitfall of the standard slow-roll calculations 
is the matching condition at "horizon crossing" k — 
0(l)aH where one conventionally just evaluates the per- 
turbations exactly at k = aH. The 0(1) ambiguity in 
the "horizon crossing" and also the possible pro- 
longed effects around horizon crossing are accounted for 
in the general slow-roll approximation. Each Fourier 
mode in the perturbation calculations is affected by the 
inflationary dynamics over the whole period it is leaving 
the horizon, not just at one instant of "horizon cross- 
ing" . Therefore, the correlations of each Fourier mode 
can be highly non-trivial compared with the cases of the 
standard slow-roll. This makes our use of principal com- 
ponent analysis rather useful because it gives us the in- 
dependent (i.e. statistically independent) modes for the 
clearer physical interpretation of the parameter estima- 
tions as discussed in Sec. IIIII 

B. Discrete Parametrization 

Scalar perturbations alone cannot constrain V itself 
(we need tensor perturbations to do it) but rather certain 
functions of V . In this paper, we consider the physical 
limitations of CMB observables without tensor pertur- 
bations for constraining 2 

fV'\ 2 V" 
Gn-3(-j -2- . (1) 

With the standard slow-roll approximation, the spec- 
tral index is directly related to this quantity G n because 
n - 1 = -3{V'/V) 2 + 2V"/V. In our analysis, we use 
the general slow-roll approximation, in which case we do 
not have such a simple relation. G n is, however, still a 
fundamental quantity representing an inflaton potential 
even for the general slow-roll cases and indeed can be in- 
terpreted as the source function for n — 1 (see Appendix 
for the derivation) 

dhxV 



J^j[-k^W'(kO][-G n }. (2) 



2 In this paper, we consider only single field inflation models. Our 
analysis can, however, be extended straightforwarly to the mul- 
tiple field inflation models because the general slow-roll formula 
for multiple-component inflation has an analogous form as that 
for the single component case I Hi . 
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where £, minus the conformal time, is an integral over 
the scale factor a up to the end of inflation: 



dt' 



(3) 



W(x) is given as 



, , 3sin(2x) 3cos(2x) 3sin(2x) , s 

W ( x ) = 2x 3 15"^ ( 4 ) 



2.7: 



so that the window function —xW'(x) in Eq. i|22|) is 

3(a; 2 - 3) 



- xW'(x) = - 3(5X ' 3) sin(2x) 4 
2a; 13 

and the normalization is such that 



r°° dx 

/ ™[-xW'(x)] = l. 
Jo x 



cos(2x),(5) 



(6) 



For the case of standard slow-roll, G n can be taken out 
of the integral and integration of the window function 
becomes unity to lead to the standard slow-roll result 
n- 1 = -G n . 

We choose the fiducial model to be a flat inflaton 
potential which leads to G„ — and study how pre- 
cisely we can constrain the deviations from G„ — 0. We 
parametrize these deviations as 



(7) 



where {pi} form a discrete set of parameters and Bj(m£) 
are defined as 



1 if In & < In f < In 
otherwise. 



(8) 



CMB observables cover less than 10 e-folds in practice 
corresponding to a small portion of a whole inflaton po- 
tential. We therefore, in our analysis, parametrize the 
inflaton potential across a range greater than 10 e-folds 
£min < £ < £max to ensure it covers the horizon cross- 
ing for observable modes. We choose a discrete repre- 
sentation rather than a continuous one through smooth 
functions because, in addition to its simplicity, the dis- 
continuities in G ra (lnf;) will not show up in V due to 
the window W(k£). In practice, we cover a large enough 
range £ OT j n < £ < Cmax and choose a sufficiently fine 
discretization that the principal components of Sec. IIII 
have converged. 

With this parameterization, the power spectrum reads 
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jW(kt)Y,PiBiQnt), (9) 




k(Mpc 1 ) 



FIG. 1: The response function d\n'P(k)/dpi as the function 
of fc(Mpc' 1 ). The plots are given for four values of pi whose 
corresponding £ values are indicated. 



where 1//q is integration constant which corresponds to 
the power spectrum amplitude of the fiducial model 



P\nd = 4 = (5-07 x 10 
Jo 
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(10) 



Here "fid" denotes that it is to be evaluated for the fidu- 
cial model (i.e. pi — 0). 

In Fig. n we pl°f t ne infinitesimal response in the 
primordial power spectra to the inflation parameter 
dlnV (k)/dpi evaluated for our fiducial model. Note that 
as might be expected the response takes the form of a 
localized change in the amount of low vs high k power 
or 'tilt' surrounding horizon crossing fc£ ~ 1. 

The CMB angular power spectra are given by integrals 
over V 



2tt 



xx' 



= J d\nkT l x (k)T l x ' 



(k)V(k), (11) 



where T x (k) are the CMB transfer functions with X rep- 
resenting the CMB temperature or E polarization (i.e. 
X is T or E). For the numerical analysis in Sec. IIIII 
we shall oversample the transfer function with 3750 log- 
arithmically spaced k modes from k — 10~ 435 — 10 -135 
Mpc" 1 and 2500 k modes from k = 10~ 135 - lO" - 35 
Mpc -1 (see |5j for details). The infinitesimal response 
in the observable power spectrum to pi, which will ap- 
pear in the calculation of the Fisher matrix, is 



dC xx 



dp t 



2tt 



fid 



1(1 + 



T) J d\nkT x (k)T l 



1 

'72 

Jo 



^ max d£ 



(12) 
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We now apply Fisher matrix analysis to calculate the 
precision to which the {pi} can be constrained to see the 
physical limitations for future CMB data in reconstruct- 
ing the inflaton potential. 



III. PRINCIPAL COMPONENT ANALYSIS 

The Fisher matrix analysis is useful for forecasting the 
uncertainties in parameter estimation, and we apply it 
to extract the information on an inflaton potential from 
Ci . Assuming a Gaussian distribution for C; , the Fisher 
matrix is given by [22^ 



F„ 



^ ^ dp u 

1=1 X,Y ^ M 



Cov 



-\c?cr 



dcr 



(13) 



where X, Y represent CMB temperature, E polariza- 
tion or their cross correlation, and Cov^ 1 is the in- 
verse of the covariance matrix of the power spectrum Ci . 
{p^} consists of the inflaton potential parameters {pi} 
and the cosmological parameters which are marginal- 
ized over: the dark energy density, equation of state pa- 
rameter, reionization optical depth and the matter and 
baryon density with their corresponding fiducial values 
fl DE = 0.72, w = -1, r = 0.17, VL m h 2 = 0.145 and 
flbh 2 — 0.024 respectively. We additionally take the 
overall amplitude of the spectrum /o as a parameter. 
This marginalization can circumvent the possible degen- 
eracies between inflationary and cosmological parame- 
ters. For quantifying the physical limitations of inflaton 
potential reconstruction, we consider ideal, cosmic vari- 
ance limited observations for both temperature and E 
polarization measurements out to l max — 2000. 

The inverse Fisher matrix approximates the covari- 
ance matrix C(« F~ ^Kramer-Rao identity) and the 
marginalized 1-er error <r(pp) for a given parameter p^ 
is 



a( Pfl ) = J(C) 



(F- 



(14) 



The errors in the individual parameters pi are in general 
large and highly correlated to one another in part be- 
cause under the general slow-roll approximation V picks 
up contributions 'around' horizon crossing not 'at' hori- 
zon crossing for each mode. This hinders the interpre- 
tation of what aspects of the potential the data will in 
fact constrain. Moreover the large number of discrete 
parameters {pi} would make a likelihood search with 
actual data unfeasible. 

We therefore instead apply a principal component 
analysis. Instead of {pi\, we consider a new set of pa- 



rameters {rrii} which are linear combinations of {pi} 

m a = [Alne] 1/2 ^5 wPl . (15) 



The normalization factor of \/Aln£, which depends on 
the discretization of ln£, is chosen so that the variance 
of m a is independent of this spacing when it is integrated 
over ln£. As long as the same convention is used for the 
signal and noise, though, the normalization is arbitrary. 
The Si a are orthonormal eigenvectors of the covariance 
matrix of the inflation parameters of interest C\j (which 
is a sub-block of C, 



)J,U) 
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Cij Sja 



1 Jn 

AlnC a 



(16) 



Therefore, in our normalization the principle compo- 
nents Sia/y/ AlnC are normalized to unit variance when 
integrated over ln£ 



Si, 



VATiT? 



\/ATn~£ 



1. 



(17) 



The {m a } are orthogonal (i.e. statistically independent) 
and their covariance matrix becomes a diagonal matrix 
with the elements 



(m a m b ) 



o 2 a 8 ab 



(18) 



We consider the parametrization of G„(ln£) as given 
in Eq. for the £ range of 10" 1 < £/Mpc < 10 4 , 
and divide it into 200 equally spaced bins in ln£ with 
200 parameters {p{\ as given by Eq. Q. The principal 
component decomposition represented by Eq. I|15fl pin- 
points the directions in the parameter space of {p{\ with 
the smallest variances. 

These normalized modes Sia/s/ A ln£ are shown in 
Fig. |21 which plots the five most tightly constrained 
principal components of the covariance matrix of the 
inflationary parameters (with the cosmological param- 
eters marginalized over) as well as 10th and 15th ones. 
Fig. |3| shows the rms error a a on the mode amplitudes. 
Since the principal components are all normalized to unit 
variance their joint constraining power can be estimated 
from the cumulative variance 



1 



y- 



(19) 



Fig. |3 shows that only the first handful of modes con- 
tribute to the cumulative variance, so we expect these to 
be the only ones worth including in an analysis. In the 
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FIG. 2: First five, 10th and 15th best principal components 
as the function of the minus conformal time £. 
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0.001 




Comparing the values of £ in Fig. [21 and Fig. we 
can infer that the first eigenmode represents the local 
tilt around k ~ 0.1 Mpc -1 , which makes sense since 
the tilt or spectral index would be the best constrained 
aspect of V(k). This is what we would expect because 
G n (In £) , which we parameterize here, represents the tilt 
as we discussed in Sec. 

Indeed, if we marginalize over the tilt n as well, the 
first eigenmode of Fig. |21 disappears, as shown in Fig. 0] 
We can also infer that the second eigenmode, which has 
compensating contributions of opposite sign represents 
the local running (i.e. the difference or change in the 
local tilt). 




10 100 
% (Mpc) 

FIG. 4: The best principal components when the tilt n is 
marginalized over. To be compared with Fig. |2] where n is 
not included in the marginalized cosmological parameters. 
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FIG. 3: RMS error on each principal component as the func- 
tion of mode number. Also shown is the square root of the 
cumulative variance defined in Eq. 1191 . 



Although these statements remain true for all models, 
understanding how constraints on the principal compo- 
nents translate into more direct statements about spe- 
cific models is more difficult. We illustrate their utility 
in the next section with examples. 



next section, we will look at some specific models and 
see that, practically speaking, this holds true: all the in- 
formation that can be profitably used to compare theory 
with observations is contained in the first few modes. 

Note that £ is given in comoving Mpc and its relation- 
ship to the number of e-folds from the inflationary era 
to the present remains undetermined and heavily model 
dependent. However we can see that the first five eigen- 
modes are sensitive to the range of 10 < £/Mpc < 500 
which covers ~ In 50 ~ 4 e-folds. This order of magni- 
tude is consistent with our expectation from ~ In l ma x = 
7.6 and the high cosmic variance of the low I modes. 



IV. EXAMPLES 

We now apply the above formalism to two simple ex- 
amples. These examples serve two purposes: (i) they 
demonstrate how a generally scale-dependent n — 1 is 
mapped onto the modes defined above and (ii) they il- 
lustrate that only the top five or so modes are likely to 
have signal to noise greater than one. This means that, 
for the purposes of comparing with CMB observations 
at least, one can parameterize inflationary models with 
only five numbers, the amplitudes of the leading modes. 



G 



The first example is a trivial one: suppose n — 1 is 
a constant, independent of scale. For concreteness, we 
choose n = 0.95. Then G n = 0.05; thus all the p^s de- 
fined in Eq. Q are equal to 0.05. Although G n in this 
model is trivial, the eigenvalues m a - which are the con- 
volution of ^ with the eigenvectors depicted in Figure |2] 
- are not. Figure [5] shows the first thirty of these. 
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FIG. 5: The first 30 eigenmodes for a potential which pro- 
duces constant spectral index n = 0.95. Error bars are the 
eigenvalues a a ; note that only the first few modes have signal 
to noise greater than one. 



This example shows that every inflation model can 
be expressed as a set of m a . Because each m a is in- 
dependent, we can obtain the total signal to noise by 
calculating \/J2i m i I °f • We found, for this example of 

n - 1 = -0.05, J2Z°i m2 J^ = 1102 which S ives us thc 
estimate of n-1 » -0.05(l±l/\/Tl02) = -0.05±0.0015. 
This can be compared with the variance of n obtained 
directly from the Fisher matrix (where n is taken as 
the only inflation parameter). In this case, we find 
cr(n) = 0.0015, in perfect agreement with the more gen- 
eral approach. We also found that including only the 
first five modes in the calculation of total signal to noise 
gives J>2i=i m i/ (T i = 1088. Hence, as long as the top five 
modes are retained, we retain essentially all the informa- 
tion about n. If we knew that n were constant of course, 
it would be better to use the constant n as the variable 
to be compared with data. 

For more general cases including the case of a variable 
V", the proper way to fit the data would be to allow 
the best constrained optimal eigenvalues, more specifi- 
cally five optimal ones, to be free parameters. Our next 
example demonstrates that this does in fact work. 

Consider a perturbation at = </>o in an inflaton po- 
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10 100 
\ (Mpc) 

FIG. 6: The change of G n as the function of £ for a potential 
with a bump given by Eq. 1120 ^ for the case of A = 0.1, c = 
lQ' G ,v = 100 and Co = 30. 



tential represented by a smooth bump 

V((/>) = VoeW-^ [1 + ce"" 2 W~^) 2 ] (20) 
where |A| < 1, |c| < 1, \v\ > 1. We now have 



d\n(aH) 



V 



1 + 



V' 



-X to the leading order (21) 



so that <j> ~ A In £ (an integration factor is scaled into 
£ such that O = Aln£ -) Eqs. (|2U|) and ipZTjl can be 
substituted into Eq. I|24|) to obtain G n as a function of 
£. This is shown in Fig.HJlfor the parameters A = 0.1, c = 
10~ V = 100 and £ = 30. 

Convolving the {pi} shown in Fig. with the eigen- 
modes in Fig. |3 we obtain the {m a } for this model; these 
are shown in Fig.[7| We see that only the first few modes 
have signal to noise bigger than unity. They therefore are 
all that need to be retained to capture the information 
necessary to distinguish this non-trivial potential from 
a perfectly flat potential. This is reasonable because, in 
this example, the deviation of G„(ln£) from zero occurs 
around the range of £ to which only the well-constrained 
modes are sensitive. 



V. CONCLUSION AND DISCUSSION 

We applied a principal component analysis to the co- 
variance matrix to investigate the precision with which 
CMB observables can ultimately constrain the inflaton 
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FIG. 7: The corresponding first 30 eigenmodes for the po- 
tential with a bump given by Eq. 1201 for the case of 
A = 0.1, c = 10" 6 , v = 100 and Co = 30. 



potential, accounting for the uncertainties in cosmologi- 
cal parameters. We parameterized the inflaton potential 
using 200 parameters in the framework of the general 
slow-roll approximation without specifying any partic- 
ular inflationary scenario. We showed that considering 
five principal components covering about four e-folds of 
conformal time will be sufficient to obtain almost all of 
the information on the features of the inflaton poten- 
tial when the deviation from the flat potential is within 
|G n (ln£)| < 0.05 in the scales CMB observables are sen- 
sitive to. Wc found the first principal component repre- 
sents the local tilt and the second one the local running 
in the region around k = 0.05 — 0.1 Mpc -1 where cosmic 
variance limited CMB observations out to I = 2000 are 
most constraining. 

We, however, should keep in our minds that the well 
constrained local "tilt" and "running" around k ~ 0.1 
Mpc -1 does not determine the properties of V(k) or 
V(4>) more globally or assess the validity of either a con- 
stant tilt or constant running. This makes measuring 
V(k) at high k from non CMB sources more interest- 
ing since in the general context, the CMB says nothing 
about n(k) beyond the compact support region of our 
best eigenmodes. Our formalism is general and can be 
easily extended to multiple data sets. For example, with 
small scale power spectrum information it may be inter- 
esting to marginalize a constant running as we did the 
tilt and construct the best eigenmode for the running 
of the running which can become rather important well 
beyond the few e-folds of the CMB range in the near 
future. 

The principal component analysis helps clarify the 
particular aspects of the inflation potential that are 



well constrained. It is also a useful technique to find 
the optimal way of extracting information on the 
inflaton potential in a model independent manner in 
the analysis of actual data, say, via the Markov Chain 
Monte Carlo method once the polarization and higher 
angular resolution data becomes available in the future. 
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Appendix: General Slow-roll Approximation 

Under the general slow-roll approximation, the power 
spectrum for a single field inflation model is given as 



lnV(lnk) = f ^-[-HW'(kO] 
Jo £ 



1 1 ^ 2 /' 

ln 7^ + 3/ 



d4 



-kZW'(kO}G(lnO, 



(22) 



The window function is defined in Eq. ©. The depen- 
dence on the dynamics of the inflaton field is encoded 



/(InO 



H ' 



(23) 



where H is the Hubble constant and cb — d<fi/dt is the 
time derivative of the inflaton field. Here and below 
primes denote derivatives with respect to the argument, 
e.g. /' = d//dln£. £ = {aH)- x [l+0(H / H 2 )], so, to the 
leading order, / = 2tt<P/H 2 which is just an inverse of the 
familiar comoving curvature perturbation. The ampli- 
tude of /'/ / is small in the general slow-roll approxima- 
tion as well as in the standard slow-roll approximation. 
/'// can, however, vary rapidly in the general slow-roll 
approximation in contrast to the case of the standard 
slow-roll approximation. As we can see from the inte- 
gral over d In £ in Eq. (|22|l , the power spectrum receives 
contributions from a range of £ for a given mode k, not 
only the one moment of horizon crossing as in the stan- 
dard slow-roll formula. This is because a given Fourier 
mode is affected by the possibly prolonged inflationary 
dynamics while it is leaving the horizon, not just at one 
instant; these details are specified by the properties of 
the window function — fc£VF'(fc£) which is shown in Fig. 
|SJ The mode oscillates rapidly inside the horizon (when 
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FIG. 8: The window function for the general slow-roll for- 
mula — k£yV' (k^) as the function of — ln(fc^). Large fc£ corre- 
sponds to earlier times, when the mode of interest is within 
the horizon, so time flows from left to right. 



k^ » 1), and one would need a sudden event to affect 
it. The window function starts to vanish for k^ < 1 
once the mode leaves the horizon because it freezes out 
and becomes hard to affect. For a slowly varying /' as 
in the case of the standard slow-roll, the window func- 
tion is not resolved and just appears as a delta function. 
Then, Eq. © reduces to V = l// 2 = (H 2 /2n<p) 2 , the 
standard result G n in the body of the paper is 

nothing but the derivative of the integrand in the general 
slow-roll formula Eq. I|22[l as follows 3 



G(ln0 




to the leading order 



(24) 



which is true even if V" changes with £ Let us 

now derive Eq. 1(5} . Covering a large enough £ range 
< £ < imax to ensure it covers the horizon crossing 
for observable modes of our interest, the right hand side 
of Eq. 1)22(1 can be reduced to 

\nV = J"~ ^ [-kZW'ikQ] G(ln0, (25) 



FIG. 9: The window function W(k( i ) as the function of 



with negligible contributions from £ outside the limits 
of integration which can be justified from W'(fc£) — * 
for the large enough £ max and the small enough £ m i„. 
Integrating by parts, 



\nV = 



-[W(kOGQnO] U - 



(26) 



W(k£) as the function of — ln(/c£) is shown in Fig. [5J 
W — > for the large enough £ max and W — + 1 for the 
small enough £ m i„ so that 

lnP = G(m£ min ) • 



where the integration constant G(ln£ m i n ) gives the nor- 
malization amplitude of V . 



3 We are here assuming a smoothly varying V(4>), see also 1171 . 
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